Bulk viscosity in kaon-condensed color-flavor locked quark matter 



Mark G. Alford,i'B Matt Braby,i'[l| and Andreas Schmitti- 

^Department of Physics, Washington University St Louis, MO, 63130, USA 
^Institut fiir Theoretische Physik, Technische Universitdt Wien, IO4O Vienna, Austria 

(Dated: September 3, 2008) 

Color-flavor locked (CFL) quark matter at high densities is a color superconductor, which sponta- 
neously breaks baryon number and chiral symmetry. Its low-energy thermodynamic and transport 
properties are therefore dominated by the H (superfluid) boson, and the octet of pseudoscalar 
pseudo-Goldstone bosons of which the neutral kaon is the lightest. We study the CFL-if" phase, in 
which the stress induced by the strange quark mass causes the kaons to condense, and there is an 
additional ultra-light "74""" Goldstone boson arising from the spontaneous breaking of isospin. We 
compute the bulk viscosity of matter in the CFL- if" phase, which arises from the beta-equilibration 
00 ' processes if" H + H and if" + H ^ H. We find that the bulk viscosity varies as T'^ , unlike the 

, CFL phase where it is exponentially Boltzmann-suppressed by the kaon's energy gap. However, in 

' the temperature range of relevance for r-mode damping in compact stars, the bulk viscosity in the 

' CFL-A'° phase turns out to be even smaller than in the uncondensed CFL phase, which already has 

a bulk viscosity much smaller than all other known color-superconducting quark phases. 

^ ■ PACS numbers: 12.38.Mh,24.85.-|-p,26.60.-|-c 
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INTRODUCTION 



The color- flavor locked (CFL) phase of quark matter is the densest predicted state of matter [ij; it may occur in 
nature, in the core of compact stars, which are expected to reach several times nuclear saturation density. Quark 
matter is described by the theory of Quantum Chromo dynamics (QCD), which, because of asymptotic freedom [1, 0], 
becomes weakly coupled and hence pcrturbativcly tractable at asymptotically high densities. In that regime, the 
CFL phase can be shown to be the ground state. However, at compact-star densities the coupling is strong, so first- 
principles calculations arc not possible. We therefore follow a more phenomenological approach. We start with the 
hypothesis that color-flavor-locked quark matter occurs in compact stars, and calculate phenomcnologically relevant 
, properties which could allow us to find astrophysical signatures of its presence. We will calculate the bulk viscosity, 

■ which is relevant for the damping of pulsations and (indirectly) the spin down of the star. 

', If the interior of the star is a perfect (dissipationless) fluid, then a certain class of oscillating modes, r-modes, are 

■ unstable with respect to the emission of gravitational waves [1, [Ei ■ This emission acts as a brake on the rotation 
QQ ' of the star. Since we know from observations that there arc compact stars with very large rotation frequencies, we 

i conclude that the instability must be damped. One damping mechanism is a nonzero viscosity of the fluid. Both 
^ ' shear and bulk viscosity can affect the r-modes and typically act in different temperature regimes. Therefore, it is of 
• ^ , physical interest to compute shear and bulk viscosity of various candidatephases in a compact star as a function of 
X ■ temperature. Several calculations exist in the literature, for nuclear 0, i i M M M Sand Im) eronic as 
^ I well as for unpaired quark matter [H, [13, and various color-superconducting phases |19l. [20l . [2ll . [2^ , for a review 
- - ' see [l^l- Bulk viscosity due to thermal kaons in the CFL phase has been computed in Ref. [241^ Here, we extend this 
calculation to the case of condensed kaons. 

At asymptotically large quark chemical potentials fiq, three-flavor quark matter is in the CFL phase [1^. This 
phase breaks the local color gauge group SU{3)c, giving rise to Meissner masses for all eight gluons. It also breaks the 
global chiral symmetry group SU{i)L x SU{3)h, leaving a residual global SU{3)c+l+r which contains simultaneous 
rotations in color and flavor space, hence "locking" color with flavor. Moreover, the global baryon number conservation 
symmetry U{1)b is broken, rendering the CFL phase a superfluid. For massless quarks there are thus 8-1-1 massless 
Goldstone bosons (and an additional Goldstone boson for the U{1)a symmetry, which, however, is expected to be 
explicitly broken at moderate chemical potentials). At moderate densities the mass of the strange quark nis cannot 
be neglected and the bosons of the octet associated with chiral symmetry breaking acquire masses. These masses are 
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small compared to the energy gaps of the fermions. and therefore the low-energy properties of the CFL phase can be 
described within an effective theory for the Goldstone bosons [26l . [27l [28j . We shall make use of this theory in this 
paper. If m^/iiq is large enough, it is expected that the lightest pseudo-Goldstone bosons, namely the neutral kaons, 
condense. The CFL phase with condensation is called the "CFL-X^" phase. Goldstone bosons in CFL and their 
condensation have also been studied in a different approach, using a Nambu-Jona-Lasinio model [29l. [sol. [Sll. [3^. Issl. Is^ . 

Besides giving rise to masses for the meson octet, a nonzero strange mass induces a mismatch in the Fermi momenta 
of the quarks that form Cooper pairs in the CFL phase. In fact, the strange mass induces a mismatch in any possible 
spin-zero color-superconducting phase (35| . This means that at lower densities, the particularly symmetric CFL phase 
may be replaced by a less symmetric pairing pattern. It is currently not known whether, going down in density, the CFL 
phase is superseded by nuclear matter or by a different, more exotic, color-superconducting phase. Candidate color- 
superconducting phases have Cooper pairs with nonzero angular momentum (sgI . [stI . [ssl [33 | or nonzero momentum 
pol. I4lll42|. In this paper, we shall only consider the CFL and CFL-K^ phases. 

The paper is organized as follows. In Sec. [II] we give a brief overview over the properties of the CFL-K'^ phase, 
in particular we present the low-energy excitations at finite temperature. As an application, we discuss the resulting 
specific heat of the system in Sec. Ill El The calculation of the bulk viscosity is presented in Sec. IIIII We define the 
bulk viscosity in Sec. IIII M In Sees. IIII B] and ITlI CI we collect the ingredients needed for the bulk viscosity, namely 
the kaon density and susceptibility and the rate of the processes X° ^ H + H and /v" + H ^ H. We put these 
ingredients together in Sec. IIII Dl to obtain the result for the bulk viscosity and give our conclusions in Sec. IIVI 



II. LOW-ENERGY MODES IN THE CFL-A" PHASE 



In this section, we briefly summarize the theoretical description and physical properties of the Goldstone bosons in 
the CFL-Ar° phase. More details can be found in the references given in the text. 



A. Chiral Lagrangian 

We denote the meson nonet in the CFL phase, associated with chiral symmetry breaking, by 

S = e*''/^' , (1) 
where 9 is an element of the Lie algebra of U (3). The physics of the mesons is described by the Lagrangian [2^, [23| 

C = ^Tr[(aoS + i[Ao,S])(9oSt-z[Ao,S]t)-v2a,I]a,S]t] +^detA^Tr[M-i(S + St)]. (2) 
There is an effective chemical potential given by the "gauge field" 

Ao = (iqQ - - — , (3) 

where Q = diag(2/3, — 1/3, — 1/3) and AI — diag{mu,md,nis) are the electric charge and quark mass matrices in 
flavor space, and fig is the chemical potential associated with electric charge. Moreover, from matching calculations 
at asymptotically large densities we know 

_ 21-81n2/^ _ J_ 3A^ 

~ 18 V3' ^ ^ 

where A is the fermionic energy gap at zero temperature. It turns out that the neutral and charged kaons are the 
lightest mesons. They carry flavor quantum numbers K'^ ~ sd, ~ ds, K'^ ^ su, K~ ^ us. In contrast to the usual 
mesons, however, they arc composed of four quarks of the structure qqqq as opposed to qq. The zero-tempcrature 
kaon masses and effective chemical potentials are deduced from the Lagrangian and are given by 

IIK+ = f^Q+ , = , (5a) 

= amd{ms + mu) , ml^o = amu{ms + md) ■ (5b) 

We see that with nid slightly larger than to„ the neutral kaon is slightly lighter than the charged kaon. Moreover, 
electric neutrality disfavors the presence of charged kaons. Therefore, in most of the remainder of the paper we shall 
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ignore the charged kaons. In the following, for notational eonvenience, we denote the neutral kaon chemical potential 
and mass simply by ^ and m, respectively, 

/i = /i/fO , TO = TOifO . (6) 

Condensation of the neutral kaons occurs if /i > m. Because of the large uncertainty in the quark masses and in the 
dimensionless quantity a, it is not clear whether this condition is fulfilled for densities present in the interior of a 
compact star. Using the high-density expressions in Eq. @ and inserting a quark chemical potential /x^ ~ 500 MeV 
and an energy gap A ~ 30 MeV, we estimate ~ 100 MeV, a ~ 0.01. With rud ^ 7 MeV, m„ ~ 4 MeV we thus 
obtain a kaon chemical potential ^ ~ 20 MeV and a kaon mass m ~ 4 MeV. These values suggest that the kaons are 
condensed and thus that the relevant phase to consider is the CFL-K^ phase. It has been argued that for sufficiently 
large values of the parameter mj./ the CFL-if^ phase is modified to the so-called curCFL-JC" phase which is 
anisotropic and exhibits counter-propagating currents from the kaon condensate and ungapped fermions [H, [H, \^ . 
Because of the presence of ungapped fermions. the transport properties of this phase can be expected to be very 
different from the ones in the CFL-i^" phase and similar to unpaired quark matter. In this paper, we do not consider 
the possibility of a kaon current but rather focus on the isotropic GF\j-K^ phase where all fermions are gapped. 

As we shall see in Sec. IIIIl in order to compute the bulk viscosity we need the nonzero-temperature behavior of the 
kaons. in particular their masses and excitation energies. The necessary results are provided in Ref. [46j and shall be 
briefly summarized in the next two subsections. 



B. 2PI formalism for kaons 



To obtain the Lagrangian for the neutral kaons in the presence of a kaon condensate we start from the Lagrangian 
^ and set 



5] ^ gj[(0-|-¥'6)T6+V7T7]//„ ^ 



(7) 



where Tg, T^ are Gell-Mann matrices. We have omitted the other Goldstone modes, proportional to 1, Ti, . . . , Ts, Tg, 
and have introduced a vacuum expectation value for the kaon field (j) (without loss of generality in the Tg direction) 
and fluctuations (pg, <^7 with vanishing expectation value. The Lagrangian is now expanded up to fourth order in the 
fluctuations. See Ref. [40| for details and the explicit form of the Lagrangian. The tree-level potential is 



T 1 9 T 

1 — cos — sm — 
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where we have abbreviated the effective coupling constant 

_ 4^2 



6/J 



(8) 



(9) 



In principle, one can keep all orders in 0, but for simplicity we have expanded for small values of (j). This restricts our 
analysis to small condensates. In other words, our results are quantitatively reliable if ^ is only slightly larger than 
TO and become unreliable for )! 3> to. In addition to the tree-level potential, the Lagrangian contains quadratic terms 
in the fluctuations from which we read off the inverse tree-level propagator 



2i/ifco 



-2i/ifco 



2 1,2 



vik 



(10) 



Here and in the following we denote four- momenta by capital letters K = (fco,k), where the bosonic Matsubara 
frequencies are given by fcp = —^imrT with the temperature T. The Lagrangian also contains interaction terms cubic 
and quartic in the fluctuations. For small condensates we may neglect the cubic interactions. The form of the quartic 
interactions is needed for the vertex in the kaon self-energy. 

The nonzero-temperature behavior of systems with spontaneously broken symmetries can be conveniently treated 
within the two-particle irreducible (2PI) formalism [i?!. I48ll49| : for the application of this formalism in chiral models 
see for instance Refs. [13, [HI, [52|. Naive thermal corrections would lead to unphysical results for a Bose-condensed 
system, in particular to negative energies for certain temperature regimes. Therefore, we apply the more elaborate 
2PI scheme to compute the thermal masses self-consistently. We start from the effective potential 
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1. 



Feff [0, S] = U{c^) + -Trln^-1 + -Tv[S^\<i>)S - 1] + [0, S] , 



(11) 
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FIG. 1: Diagrammatic representation of tiie two-loop approximation for V2 and corresponding kaon self-energy. The lines 
represent the full kaon propagator S, to be determined self-consistently. 



where the trace is taken over momentum space and over the two-dimensional space given by the two degrees of 
freedom of the kaon (corresponding to Tg and Ty or and The effective potential is a functional of the vacuum 
expectation value (j) and the kaon propagator S. The term V2[(f>, S] is the sum over all 2PI diagrams. We employ the 
two-loop approximation for this infinite sum and only consider the "double-bubble" diagram, see Fig. [1] In this case 
V2 docs not explicitly depend on (j). 

One can now determine the condensate and the thermal kaon mass self-consistcntly. This is done by solving the 
stationarity equations 







dU 



iTr 
2 



where S = is the kaon self-energy, see Fig.[T] With the ansatz for the full inverse propagator 



the stationarity equations read 



-ki + vie + Ml 

2ifikQ 



(Af 2 _ ^2)^ 



—2ifikQ 

^ki + vie + Mi 



2a(j)'^ + 2aI{M^,T) . 



(12a) 
(12b) 

(13) 



(14a) 
(14b) 



Here, M is the temperature dependent mass, related to M± in the propagator ^ by RP = {Ml + Ml) /2. The 
other combination 5M^ = {Ml — M'^)/2 can be eliminated since the second matrix component of Eq. (jl2bp implies 
(5A/^ ~ ad)^} We have abbreviated 



I{M^,T) 



(27r)3 Ek 



with 



and the Bose distribution function 



(15) 



(16) 



(17) 



The kaon excitation energy is a pole of the propagator (|13p . Its form is crucial for the thermodynamic and transport 
properties of the GF\j-K^ phase and will be discussed explicitly in Sec. Ill C[ see Eq. p^ . The second pole of the 
propagator corresponds to the excitation. This excitation has a (temperature-dependent) energy gap which is 
larger than the excitation by at least 2/i. We will neglect the mode in our analysis, because ^ is expected to 



^ In fact wc have (5M^ = aij)^ /{I + a J) with J being a momentum integral similar to / in Eq. I I15I I. We may neglect the term a J . This 
removes a very small mass for the Goldstone boson, which arises as an artefact in the 2PI formalism at finite temperature. For details 
and the expression for J see Ref. [46| . 
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be on the order of ~ lOMeV, so the K'^ would only play a role at temperatures in the tens of MeV range. This is 
close to the critical temperature of the CFL phase itself, and is not physically relevant because compact stars cool 
below such temperatures within minutes of their formation. Our results are therefore reliable at temperatures of order 
1 MeV or lower. In some cases we will continue our expressions to higher temperatures: this represents a theoretical 
exercise in studying the properties of the pure K'-' condensate, not a physical prediction. 

For the calculations of the specific heat in Sec. Ill El as well as for the density and susceptibility in Sec. HUB] we need 
an explicit expression for the pressure P which is the negative of the effective potential P = — V^ff . In the given 
approximation we can write P as a function of ifP and A/^, 

P{4>',M') = - tI ^ In (l - e-^/^) + ^(M^ - - 2a02) /(A/^, T) - ^P{M\ T) . (18) 

The pressure in the physical state is then given by inserting the values for (p and AI at the stationary point. 

C. Kaon excitation energies 

For the calculation of transport properties we are interested in the kaon dispersion relations which are the poles 
of the pro pag ator (|13p . Dropping the i?" excitation, see remark below Eq. (|17p . the only relevant kaon dispersion 
relation is 1461 



1/2 



for T <T, 



C -I 



efc=< ' ' (19) 



^Jvlk'^ + M^{T) - ^ for T > , 

with M{T) determined from the stationarity equations (jl4p . We have eliminated </> from the expression for by 
using Eq. (|14ap . The critical temperature of the second-order phase transition for kaon condensation is denoted by 
Tc, i.e., for temperatures below (above) Tc we have ^ (</) = 0). The value of Tc is estimated to be at least 
of the order of tens of MeV [i^. This shows that for all temperatures relevant for compact stars, the kaons will 
be condensed (provided the parameters m and /i arc such that there is condensation at zero temperature). The 
self-consistent treatment ensures that for all temperatures M{T) > fj,, and thus the energies given in Eq. p9)) are 
real and positive as they should be. Moreover, from Eq. (fTO|) we see that £^=0 = for T < T^. This shows that 
there is a massless Goldstone mode associated with kaon condensation.^ It is the presence of this gapless mode that 
causes thermodynamic and transport properties of the CFL-/\ ° phase to differ from those of the CFL phase. Strictly 
speaking, this Goldstone mode also has a small energy gap. This energy gap does not show up in our treatment 
within the effective theory given by the Lagrangian ([2]). The energy gap is rather induced by the weak interactions 
and can be estimated to be in the keV range [5^ . Most of the temperature scales we show in this paper do not include 
temperatures below this range. Therefore, we shall neglect this effect and speak of an exact Goldstone mode with 
massless dispersion relation given by Eq. (|19p . 



It is instructive to expand the kaon dispersion below Tc for small momenta. 



<fc + r..2.^", ..2i3 <fc - (20) 



For small momenta, the dispersion is linear in k. However, for T Tc, the coefficient in front of the linear term goes 
to zero since M{T Tc) — > ^. Consequently, around Tc the quadratic part of the dispersion becomes important. 



D. One- loop approximation for the superfluid mode H 

Next we turn to the H boson which is associated with the spontaneous breaking of U{1)b- This Goldstone boson 
can be described by the effective Lagrangian 

^ = \ido^)' - ^(v^)^ - ^ao^a^^av + j^4(9,^9^vr , (21) 



^ We shall in the following use the term K'^ for this Goldstone mode in the CFL-ii'" phase as well as for the member of the meson octet 
arising from chiral symmetry breaking in the CFL phase. This is the most convenient terminology but one should keep in mind that 
these two excitations have very different dispersion relations. 
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n = 



FIG. 2: Diagrammatic representation of the one-loop self-energy 11 for the superfluid mode H . The lines represent the tree-level 
H propagator Do- The vertex, here denoted by a square, is given in Eq. ((24}. 



where vh = l/'v/S- For notational convenience, we shall from now on abbreviate both vh and by v, 

v = VH = v-^ . (22) 
For the calculation of the bulk viscosity we need the one-loop self-energy of the i?, given by the diagram in Fig. [21 

n(^) = W^T7ll ^(^' K)D^{K)D,{P K) , (23) 

where the vertex (squared) is 

F{P, K) = [pai2K ■P-K^) + ka{P^ - 2K ■ P)f , (24) 
and the inverse tree-level propagator is Dq^{K) = fcg — v^k^. We need the imaginary part of the retarded self-energy 



Imn(P) = ^77^ J2 I t A U F(.eivk)f{e,vk)f{e2v\p - k|)^(po - e,vk - e,v\p - k|) , (25) 

8lMg /(Po) 7 (27r)3 4v2fc|p- k| 



where F{eivk) = F{P, K)\ko=eivk- The angular integration can be performed exactly. One finds that the self-energy 
assumes different forms depending on the sign of pq — vp, 

imn(p) = n+(p)e(po - vp) + n-{P)e{vp-po) , (26) 

where 

n^(^) ^ ^^5rrT-77;rT / " dkg'{P,k)f{vk)fipo ~vk) , (27a) 
n-(P) ^ -^^^ / dkg'-{P,k)f{vk)[l + f{vk~po)]. (27b) 

Here we have abbreviated 

5(P, k)=pl~ - 3(1 - v'^)vk{po - vk) . (28) 
Denoting the one-loop H propagator by 

D-\P)=D^\P) + Tl{P), (29) 
and neglecting the real part of 11 we have an approximate form of the imaginary part of the propagator, 

^"-^(^)- (pg^.V%mn2(P) ' 

which we will need in the calculation of the bulk viscosity in Sec. IIIIl 
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E. Specific heat 

The self-consistent formalism from Ref. j46| . summarized above, provides us with the tools to compute thermody- 
namical quantities of the CFL-i^T^ phase. As a first application, we shall discuss the calculation of the specific heat. 
This is of physical relevance for instance for the cooling behavior of compact stars; see for instance Refs. P, [s^, [s^l 
for the specific heat of other quark matter phases. We will not need the result for the calculation of the bulk vis- 
cosity. However, the calculation shows in an exemplary way how to compute thermodynamical quantities in the 2PI 
formalism and will then in Scc. lIIIBl bc applied to the calculation of the kaon susceptibility. 

The definition of the specific heat at constant volume cy is 

cv = r-, s^ — , (31) 

where s is the entropy density. It seems straightforward to take the second derivative of the pressure in Eq. with 
respect to the temperature. However, the self-consistent treatment complicates this procedure. Since the pressure is 
a function of the self-consistent quantities (j> and M, the implicit dependence of (j) and M on the temperature and the 
constraint of the stationarity equations have to be taken into account. We present the details of the calculation in 
Appendix The result is 

, _ 1 /• el 1 T dl dA-P 

^ 2 J (27r)3 r2 cosh f - 1 2 (9T 9r ' ^ ' 

where / is the integral defined in Eq. ([15]), and dM^/dT is given in Eq. (jA8bp . From this expression one obtains the 
discontinuity Acy of the specific heat at the critical point. 

For all interesting parameters, the first term on the right-hand side of Eq. (j32p is dominant and Acy is small, as we 
shall see from the numerical results. For small temperatures only small momenta of the kaons contribute. Therefore 
we can use the small- momentum kaon dispersion from Eq. (|20p to approximate the specific heat at small temperatures. 
The linear part of the dispersion produces a contribution cubic in the temperature, 

15^3 Y 



where we have approximated the thermal mass by its zero-temperature value, ~ 2fi^ — m? . The dispersion 
(I20[) shows that at the critical temperature the linear term vanishes and the dispersion becomes quadratic in the 
momentum. In this case, the contribution to cy goes like T^/^, 



quad 



aT3/2 (3Ai2-m2)9/4 



211/4^2^3 ^3 ' COShx^ 



oo 



dx — -~ 4.46. (35) 



The kaon contribution has to be compared to the H contribution. Approximating the dispersion by the linear behavior 
ufc, cf. the inverse tree-level propagator below Eq. (j24p . we find for low temperatures 

By comparing this expression with Eq. p4p we conclude that the low-temperature contribution of the has the 
same T dependence as the H contribution, but that the contribution is always larger by a numerical factor which 
depends on the zero-temperature kaon mass and chemical potential. We show the full numerical result of the kaon 
specific heat and its comparison with Cy", Cy"*^, and Cy in the left panel of Fig. [3l We see that there is a significant 
temperature regime in which Cy^'^ is a good approximation to the specific heat, i.e., cy behaves as T^/^ as opposed 
to . The presence of this temperature regime depends on the parameters. Here we have chosen the kaon mass and 
chemical potential such that this regime is relatively large (/i only slightly larger than m). For larger values of /i, this 
regime is typically smaller or hardly visible (see right panel of Fig. [3]) . Note also that our temperature scale extends 
to very small temperatures in the keV regime. In this regime we expect the small mass of the kaon, which is not 
present in our formalism, to play a role, see remark below Eq. (jl9p . 
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TABLE I: Collection of thermodynamic properties of neutral kaons in the CFL-_fL° and CFL phases. Here "high T" means in 
particular T > Tc, i.e., there is no CFL-A''^ phase at these temperatures. The zero-temperature kaon mass is denoted by m, 
and (5m = m — n with the kaon chemical potential fi. In the CFL-_R'° phase we have 6m < 0, and we show the low-T results to 
lowest order in Sm/m which corresponds to small condensates ^//ir ^ 1- In the CFL phase, Sm > 0. 



The discontinuity at the critical point is small and not visible in the plot. At large temperatures, the specific heat 
again goes like T^, however with a different prefactor than at low temperatures. In Table |T] we have listed the low-T 
and high-T behavior for the kaon specific heat together with the density and susceptibility, to be computed in Sec. 
ImBl In the table we also show the behavior for the case of no kaon condensation, i.e., for m > ^. 

So far we have ignored the charged kaons. Due to the electric charge neutrality condition their number density 
is suppressed at small temperatures (46j . However, in the isospin symmetric case they give rise to an additional 
Goldstonc mode [i^, [H, [11] which in principle gives a large contribution to the specific heat. To illustrate this 
contribution we ignore the neutrality constraint for simplicity and compute the specific heat for a fixed charged kaon 
chemical potential fxj^+ . For low momenta and exact isospin symmetry fj, = j = the additional Goldstone 
mode has a quadratic dispersion, = u^fc^/(2/i^). Consequently, its contribution to the specific heat goes like T'^/^, 
just as the contribution from Eq. (j35p . 

We may use the formalism of Ref. [46t to compute the thermal masses in the two-component system of neutral and 
charged kaons and derive the specific heat analogously to the above one-component system. The numerical result is 
shown in the right panel of Fig. [31 For this plot we have chosen the kaon masses and chemical potentials such that 
the isospin symmetry is almost exact. We see that there is a large contribution from the charged kaon mode which 
disappears for sufficiently small temperatures. This is due to the small isospin symmetry violation in the parameters 
and hence a small charged kaon energy gap. 



III. BULK VISCOSITY IN THE CFL-i^" PHASE 



We now turn to the main part of the paper where we use the results of the previous sections to compute the bulk 
viscosity in the CFL-A'" phase. 



A. Definition of bulk viscosity 

The definition of the bulk viscosity and the derivation of its expression in terms of the kaon rate can be found in 
Ref. [2^ . Therefore, here we only briefly summarize the most important relations and the underlying physics. In 
general, a relativistic superfluid may have more than one bulk viscosity, related to stresses in the superfluid flow with 
respect to the normal flow [ll|, [gO, [Mi- Here we neglect this effect and compute the bulk viscosity related to the 
normal fluid. 

We arc interested in a system with volume Vq which undergoes a volume oscillation with amplitude dVo ^ Vb and 
frequency uj, 

V{t) = Vn+5Vo cos ujt. (38) 

In the astrophysical setting, the oscillations are local volume oscillations with an (inverse) time scale typically of the 
order of the rotation frequency of the star which can be as large as w/ (27r) ~ 1 ms~^. The periodic change in volume 
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FIG. 3: (Color online) Various contributions to the specific heat as a function of temperature in the CFL-K'^ phase. In 
both panels we show the specific heat of the if" mesons as a solid black line, and the low-temperature specific heat of the H 
mesons as a dashed-dotted (blue) line. These coincide at large T because the low-T approximation for the H is identical to 
the large-T approximation for the K'', see Eq. (|36|l and Table [T] In the left panel we show low-temperature approximations 
for the contributions of the linear (steeper dashed line, Eq. (|34p ) and quadratic (shallower dashed line, Eq. (|35pl parts of the 
kaon dispersion. These are good approximations in different temperature ranges. In the right panel we show the contribution 
of the mode for the case of an almost exact isospin symmetry (curved solid (red) line) and the approximation from 
Eq. (|37p (straight dashed (red) line). Parameter values: for the left panel we have used a kaon mass ?n = 4MeV, and a kaon 
chemical potential fj, = 4.01 MeV; for the right panel we have used m — 4MeV, ^ = 4.501 MeV for the neutral kaons and and 
— 4 MeV, — 4.5 MeV for the charged kaons. These parameters lead to critical temperatures of Tc ~ 8 MeV (left) 
and Tc ~ 39 MeV (right). 



induces a change in density. This change in density, in turn, may induce a change in the chemical composition of 
the matter. In CFL-if° matter, there will be a change in the strangeness content, described by an induced nonzero 
(5/i, i.e., the equilibrium kaon chemical potential /i is shifted to the nonequilibrium value -I- dfi. The response of 
the system is to create (or annihilate) kaons in order to reequilibrate. The dominant processes considered here that 
change kaon number (and thus strangeness) are 

^ H + H , (39a) 
K° + H ^ H . (39b) 

These processes have also been considered in Ref. [IJ] for the case of uncondensed kaons. If the kaons are condensed 
there is also a cubic interaction, induced by attaching one leg of the quartic interaction to the condensate, which 
directly moves strangeness between the condensate and the thermal kaon gas. Here we neglect these processes since 
they are proportional to the condensate. This is consistent with our 2PI treatment which is only valid for small 
condensates, i.e. (/x — m)/m <C 1. 

The processes ([5^ induce, in response to the external oscillation V{t), an oscillation in the kaon chemical potential 
fj-{t). If the external oscillation and the system's response are out of phase, there will be dissipation, resulting in a 
nonzero value of the bulk viscosity. The bulk viscosity is maximized if the external oscillation and the rate of kaon 
production (annihilation) is on the same time scale. In this sense, bulk viscosity is a resonance phenomenon. In fact, 
it is the exact analogue of an electric circuit with alternating voltage that responds by an induced alternating current 



The definition of the bulk viscosity is 



where the dissipated power is given by 



._2_KL^ (40) 



^^ = -l^.Im^P. (41) 

dt 2Vq ' 

The imaginary part of the complex amplitude SP is given by 

ImdP = nlmSfi + nqlmSfiq , (42) 
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where the complex amplitudes 5fi, Sfiq account for the oscillating change in the chemical potentials, 

fi{t) = ^i + Rc{6fie''^*), (43a) 
fig{t) = ^ig + Re{5nqe'^*). (43b) 

Here, Uq and /Zg are the quark number density and chemical potential, while n and ^ are the corresponding kaon 
quantities. Note that the effective potential ([TT|) also depends on the quark chemical potential fiq. In particular, 
the quark density is nonvanishing and there is an induced oscillation also in fiq. This effect, however, is small. In 
Appendix [B] we present the general derivation of the bulk viscosity, taking into account the quark density effect. Here 
we proceed by only keeping the kaon terms. In this case, there is a single differential equation for 

dn dii dn dV ^ , , , , 

The left-hand side of this equation simply means that the kaon density is a function of the kaon chemical potential, 
hence the time dependence of n can be expressed in terms of the time dependence of /i. The right-hand side expresses 
this change in terms of the volume change (first term) and the kaon rate (second term) . The kaon rate F^^ o is defined 
as the change in kaon number per time and volume due to the processes p9|) . For the microscopic definition of F^^^o 
see Eq. (|58p in Sec. IIII CI For sufficiently small 8^ we can approximate F^o by 

Fj<-o ~ -ARe((5/ie*"*). (45) 

This defines the factor A to be positive (for 5^> Q, the system responds by annihilating kaons, hence F^^-o < 0; on the 
other hand, if < 0, the system responds by creating kaons and thus F^fo > 0; in both cases we thus have A > 0). 
Inserting Eqs. ([38]), (|43ap . and (gSl) into Eq. ([44]) and inserting the resulting solution for ImS/j, into Eqs. (|42]), ((4T|). 
and (|40p yields the bulk viscosity 

C = 4 2^a/ (46) 



where we have denoted the kaon number susceptibility 



dn 



In the following we shall evaluate the bulk viscosity ()46p . To this end, we shall compute the equilibrium density n 
and susceptibility x in Sec. IIII Bl the kaon rate A in Sec. IIII C[ and put the results together in Sec. IIII Dl 



B. Kaon density and kaon susceptibility 



The kaon density is given by the negative of the derivative of the effective potential pT|) with respect to the chemical 
potential. 



dU_ 



-Tr 



dix 



(48) 



Here we have only taken the explicit derivative with respect to jjL. Implicit dependencies through the self-consistent 
condensate (j) the kaon mass M drop out of the density when we take the value of n at the stationary point. 
Inserting the tree-level potential ^ and the propagators ([TU)) and ([T^ yields (see also Ref. [i^) 



1 



1 da 



(27r)3 



(49) 



with the integral / defined in Eq. psp and the kaon excitation energy from Eq. (|19p . At zero temperature, the 
integral over the Bose distribution function (as well as the integral /) vanishes, and the density is nonzero only 
because of a nonzero value of 0; in other words, all kaons are condensed. Note that the effective coupling a, see Eq. 
©I depends on the kaon chemical potential. With the zero-temperature value of the condensate (j)^ = (/.j^ — m^)/a 
we easily obtain the low-temperature limit of the density. 



2) {2^l^'+m^) 



(4^2 



(50) 
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FIG. 4: Kaon density (left panel) and kaon susceptibility (right panel) in the CFL-T^" phase as a function of temperature 
(solid lines: full numerical result; dashed lines: low-T and high-T approximations from Table |T]|. Both quantities have nonzero 
values at T = 0. We have used a kaon mass m = 4MeV, and a kaon chemical potential ^ = 4.3 MeV. This leads to a critical 
temperature Tc — 49 MeV. At T = Tc, the density has a kink (not visible on the plotted scale), and the susceptibility is 
discontinuous. 



The expansion of this limit for small condensates and the opposite limit for high temperatures are shown Table ID 
The two limits are compared with the full numerical result in the left panel of Fig. 

The derivation of the susceptibility x is more complicated. The susceptibility is the second derivative of the pressure 
with respect to the chemical potential, and the dependence of and M on fi does not drop out. Hence we have to 
compute X analogously to the specific heat, which is also a second derivative of the pressure. Details of the calculation 
are presented in Appendix [Ul The result is 



Xo + Xi + 



1 

2T 



1 



Xl + 7^ 



1 

2T 



(27r)3 cosh ^ 



1 



(27r)3 cosh^ 



1 



for T<Tc 



for T>Tr 



(51) 



where xo^ Xij ^'iid Xi given in Eqs. (|C4|) and (|C5 
approximating the condensate by its zero-temperature value, 

-6 u15jj^4^2 



The zero-temperature limit is obtained from xo upon 
2 ^ _ jYi'^y^^^ 



(4^2 _ ^2-) 



2^3 



(52) 



The high-T behavior is dominated by the momentmir integral in Eq. (|5ip . Table U shows both limit cases. The full 
numerical result for the susceptibility and its comparison to the limit cases is shown in the right panel of Fig. [51 In 
contrast to the density and the specific heat we observe a strong discontinuity at the critical point. 



Rates of the processes A'" ^ H + H and K° + H ^ H 



The final ingredient to the bulk viscosity (jlS]) is the rate F^o due to the processes To this end, we need the 

imaginary part of the kaon self-energy S given by the diagram in Fig. O This is simply the imaginary part of the 
one-loop H propagator ((30)) multiphed by the square of the K^H vertex GdsfTrfHiPo ~ ^dsP^) 01, where 



Gds = —V2VudVusGF , 



Vds 



1 

7!' 



2 _ 



(53) 



with the Fermi coupling Gp and the entries of the CKM matrix Vud, Vus- 
be understood as the weak process d + s ^ u + u, for details see Ref. [H 



(Microscopically, the K''^H interaction can 
.) Consequently, we obtain 



ImS(P) 



GlPjfjipl-vj^p'rimlliP) 



Imn2(P) 



(54) 
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FIG. 5: Kaon self-energy E needed to compute the rate of the processes ^ H + H and K'^ + H ^ H. Thin dashed hnes 
are tree- level H propagators Do, thick dashed lines are full H propagators D, and solid lines are K° propagators. The if" H 
vertex, here denoted by a cross, is GdsfirfHipo — ""dsP^)- The H vertex, denoted by a square, is given in Eq. (|24|l . 



We can now compute the kaon production rate F^o . In the closed-time path formalism (for calculations in a similar 
context using this formalism sec for instance Refs. [13, [13, El, ) the rate can be written as 



Tifo = i I |]^[G<(P)E>(F) - G>(P)S<(P)] . (55) 



Here, the self-energies are given by 



and the propagators are 



E>(P) = -2z[l + /(po)]ImE(P), (56a) 
E<(P) = -2*/(po)ImS(P), (56b) 



G>(P) = _^[i + /(e;)]^(p„_e;-<5A.), (57a) 

G<{P) = - — f{e'^)6{p,-e'^-5^^). (57b) 

In these propagators we have introduced chemical noncquilibrium by a nonzero Sfi which enters the kaon excitation 
energy e'^ = ep{ii fj, + 6fi). In ep the self-consistent mass M is determined according to the modified chemical 
potential. This is necessary since the value of the kaon mass M is governed by the strong interaction and thus the 
mass adjusts itself instantaneously compared to the equilibration time of the weak processes p9p . The energy -I- (5/x 
in the (5-functions of the propagators can be understood from the case without kaon condensation. In this case, we 
have ep = Ep — ^, and Sfi simply cancels, + = ep (ignoring the change in the self-consistent mass). In the case 
with condensation, this cancellation is only partial. Inserting Eqs. ([56|) and ([57|) into Eq. ([55]) yields 

r^o = GL/l/i / (^^^(Po - 4 - ^A^) {[1 + /(4)]/bo) - /(4)[1 + fiPo)]} Jf„;vf ■ (58) 

From this expression it is clear that J/i = corresponds to chemical equilibrium because in this case forward and 
backward rate cancel each other, and F^^o — 0. For small 6fi we may approximate 

This shows that the sign of the rate is determined by the sign of 6fi. For positive 5fi the system reacts with annihilating 
kaons, and F^o < 0, while for negative Sfi kaons arc created, and Fj^o > 0. With the approximation ([55]) we see that, 
to lowest order in Sfi, the rate is 

, . _^ , sM/i f ^ ,,^), , ,,,,, (4 . 

Sfi T J 2ep(27r)-^ (e^ - u^p^)^ + Imn^(ep, p) 



in agreement with Ref. [2J]. We shall use Eq. ((60l) for the numerical evaluation of the rate. Before we turn to the 
results let us first rewrite the rate into a more instructive form. Inserting the H self-energy from Eqs. (j26p and (j27p 
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into Eq. (|58p . wc obtain 



{e(e; + ^^^- vp)[l + fie'p)] /, dk g\e'^ + k)f{vk)f{e'p + d> - vk) 



2Q{vp - e; - Sfi) [1 + /(e;)] / , ^^^^^^^ dfc .g^ (e^ + Sfi, p, k)f{vk) [1 + f{vk-e'^- 5fi)] 



- 6(6; + 5pi- vp)f{e;) I, &52(e; + Sti,p, k)[l + f{vk)][l + /(e^ + Sfi~ vk)] 

- 2Q{vp - e; - /, ^^^^^^ dfcg2(e; + <5/.,p, k)[l + /(«fc)]/(«fc - - <5^)} . (61) 

This expression is useful to discuss the separate processes that contribute to the rate. From the structure of the Bose 
distribution functions (interpreting / as an ingoing and 1 + / as an outgoing particle) we see that the kaon rate is 
composed of four contributions: the processes H + H ^ K'^ and H K'^ + H (first and second term, respectively) 
and the inverse processes ^ H + H and H + K'^ H (third and fourth term, respectively). As expected, we 
see that the kaon-creating processes yield a positive contribution to F^o while the kaon-annihilating processes yield 
a negative contribution. Furthermore, wc conclude that in the CFL-/v'^ phase (and for infinitesimal 5^) only the 
processes H ^ + H are allowed. This can be seen from the dispersion relation in Eq. p9p which implies < vp 
for T < Tc- Then, due to the step functions in Eq. ([6T|) . the process H + H ^ K'^ is excluded. This is of course a 
direct consequence of energy conservation. 

For small temperatures we can derive an analytic expression for the rate in the CFL-X*^ phase. The derivation is 
presented in Appendix [P] The result is 



AcFL- A'o - y Q(m, M) , (62) 



where Q{m, fj.) is a complicated dimensionlcss function of the zero-temperature mass m and chemical potential /x given 
in Eq. (|D7[) . We can simplify this expression by defining 



Sm = m ~ fj, (63) 
and expanding for small values of |(5m|/m, which corresponds to small kaon condensates. This yields 



80 GIMPhT^ 



AcFL-XO — — — — T 



l + O 



\6m\ 



(64) 



The full evaluation of the rate has to be done numerically; we have checked that the expression ([S^ is in very good 
agreement with the full result up to temperatures T < 1 MeV. The additional approximation (jM]) is accurate to 
within 30% for values of \5m\/m < 0.05. 

In the left panel of Fig. [5] wc show the numerical result for A as a function of the temperature. From Eq. (I53p we 
have Gds — —7.8 • 10^^'^McV^^; moreover, we use fn — Jtt — 100 MeV. We compare the rate in the CFL-if" phase 
with the one in the CFL phase. To realize these two cases we choose the kaon chemical potential slightly above and 
slightly below the kaon mass. 

The main features of the two rates are as follows. At very small temperatures, the rate in the CFL-if ° phase is 
larger than in the CFL phase; parametrically, the former behaves as T'' while the latter is exponentially suppressed. 
At very large temperatures, the two rates are almost identical. In the intermediate temperature regime however, 
which is the regime relevant for neutron stars, the rate in the CFL phase is much larger. The reason is a larger phase 
space for the weak process, as we shall explain in detail in the following. 

The kaon width is largest for kaons of momentum p = p such that the internal H particle in the kaon self-energy 
diagram Fig. \S\ is closest to being on shell: the H and A'*' with this momentum both have the same energy, Cp = vp, 
and the denominator of Eq. (|60p is then minimized. From the kaon excitation energies in Eq. (jl9p we conclude 

^ CFL, 

(65) 

CFL-K° . 
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FIG. 6: Left panel: Kaon decay rate A = — F^o /(S/i as a function of temperature. We compare the rate in the CFL-T^" phase 
(solid line) with the rate in the CFL phase (dashed line). We have set the kaon mass to m = lOMeV. The kaon chemical 
potentials are chosen slightly above, jj, = 12MeV, and slightly below, = 8MeV, the kaon mass (i.e., Sm = — 2MeV for 
CFL-if", Sm = +2MeV for CFL). The quark chemical potential is = 400 MeV. The critical temperature is located at 
the right boundary of the temperature scale, i.e, the entire solid curve describes a kaon-condensed system. Right panel: bulk 
viscosity ^ for the same parameters and a compression frequency of uj/{2-k) = Ims^^. 



In the CFh-K'^ phase there will be much less phase space associated with the maximum in the kaon width, because 
it occurs at zero momentum, so one expects the total kaon rate to be lower. Whether the additional phase space in 
the CFL case is available and has a significant effect on the rate depends on the temperature. For sufficiently large 
temperatures, T 3> p, where states well above p = p are thermally populated, the effect of the peak at p = p is 
negligible. Moreover, for large momenta the kaon dispersion is independent of whether there is a condensate or not; 
the dispersion then is simply linear, ep ~ vp. This is the reason why the two curves in the left panel of Fig. [S] are, 
at large T, almost on top of each other (the small difference in the curves is simply due to the different values of fi, 
not because of any qualitative difference between the phases). In other words, besides the values for p given in Eq. 
(|55|). ep = vp is also satisfied asymptotically for large p (in both CFL-i^'" and CFL phases alike). This momentum 
regime becomes available at large temperatures, and thus the rates become almost identical. For smaller (but not too 
small) temperatures the peak in the integrand in the CFL phase is responsible for the large kaon rate compared to the 
CFL-i^'" phase. For sufficiently small temperatures, T <^p, states around p = p are no longer populated. Therefore, 
the rate in the CFL- ivT^ phase, which has its largest contribution from momenta close to zero, becomes larger than 
the rate in the CFL phase which becomes exponentially suppressed. 

D. Results for the bulk viscosity 

We can now put together the results from the previous subsections to compute the bulk viscosity. We insert the 
numerical results for the density n, see Eq. and left panel of Fig. 01 for the susceptibility x, see Eq. ([51]) and right 
panel of Fig. IH and for the kaon rate A, see Eq. ((60)) and left panel of Fig. [6l into the expression for the bulk viscosity 
(H51) . We shall first analyze the features of the bulk viscosity below the critical temperature and compare the result 
with the bulk viscosity in the CFL phase. We shall use dm = m — fi a.s a parameter to distinguish between the CFL 
phase {Sm > 0) and the CFL-/v° phase {6m < 0). Then, we analyze the behavior of the bulk viscosity at the critical 
point. This of academic rather than physical interest since the critical temperature is likely to exceed temperatures 
reached in a compact star. Finally, we shall discuss the parameter dependence on Sm and compare the results with 
the bulk viscosity in other quark matter phases. 

1. CFL-K° vs. CFL bulk viscosity 

We first use the parameters of the rates shown in the left panel of Fig. [S] to compute the corresponding bulk 
viscosities of the CFL and CFL-/C° phases. The result is shown in the right panel of Fig. [B] We have chosen an 
oscillation frequency uj/{2'k) = Ims^^ which is typical for a compact star. Let us first discuss the gross features of 
the result, valid for both CFL and CFL-if" phases. As explained in Sec. lIII Al if the prefactor r? is held constant, 



15 




1.0 X 10 



1.10 



FIG. 7: (Color online) Behavior of the bulk viscosity around the critical temperature, here chosen to be Tc — 37MeV (m = 
lOMeV, ^ = 10.5 MeV), for different frequencies v = aj/(27r). For > 1 ms~^ the discontinuity becomes very noticeable. 



the bulk viscosity becomes maximal if the frequency w and the rate of the processes ([5^ are on the same timescale. 
More precisely, we have to compare oj with the effective rate A/x which appears in the bulk viscosity. As a very rough 
estimate, we read off from Fig. |4] that the susceptibility is of the order of x 10"'McV^ (using ~ 100 McV). This 
means that the kaon rate has to be of the order of A ~ 10~^'^McV^ to render the effective rate A/x of the order of 
a frequency in the ms~^ regime. We see from the left panel of Fig. [Hlthat this is the case for large temperatures, 
T > 20 MeV. And indeed we observe that the bulk viscosity has a maximum in this temperature regime. Far below 
these temperatures we have lu 3> A/x, and the bulk viscosity can be approximated by C — Xn^ / {x^ lo^) ■ We now 
discuss the CFL-X^ and CFL results separately in this low temperature regime. 

In the CFL- if ° phase, the entire temperature dependence of the low-temperature bulk viscosity is given by the rate 
A because the density n and the susceptibility x tend to constant values at low T. This can also be seen by comparing 
the two curves in the left and right panels of Fig. [51 With the low-temperature expression for the rate ((M]) and the 
values for n and x from Table U we obtain the bulk viscosity in the CFL-/\" phase at small temperatures. To lowest 
order in \5m\/m, corresponding to small condensates, we find 

In the CFL phase without kaon condensation, the rate behaves very differently at small temperatures: it is 
exponentially suppressed, Acfl exp(— (^to/T). The prefactor n^/x^ contributes an additional factor of T^, 
CcFL = T^Acfl/'j-'^- In Fig- [6] we confirm that the bulk viscosity at low T in the CFL phase has a different T 
dependence than the kaon rate. 



2. Critical behavior of the bulk viscosity 



Next we turn to the behavior of the bulk viscosity at the critical temperature where kaon condensation disappears, 
and there is a transition from the CFL-if° phase to CFL. This is probably not relevant to astrophysics, for two 
reasons. Firstly, one would have to fine tune the kaon mass and chemical potential to bring this critical temperature 
down to typical compact star temperatures. Secondly, the critical temperature for kaon condensation is of the same 
order as the critical temperature for the underlying CFL condensate itself p6| , so there may really be a transition to 
unpaired quark matter. However, from the theoretical point of view it might be interesting to compare the critical 
behavior of the bulk viscosity in the present context with the critical behavior in related systems. For example, the 
bulk viscosity in a hot quark-gluon plasma has recently been studied (64,, JiS,, 6^ , and simple models not unlike the 
one we use for kaons in this paper have been employed to study the behavior at the critical point [gJ , see also [6^ . 

With the help of Eq. (|46l) we can make qualitative predictions of the behavior of the bulk viscosity at the critical 
point. For small values of the frequency, u) ^ A/x, the viscosity behaves like C — /\- Both the density n and the 
rate A are continuous at the phase transition, hence the bulk viscosity is continuous too. On the other hand, for large 
frequencies, u) ^ A/x, the bulk viscosity behaves like C — 7t.^A/(x^<^^). From Fig. 2] we know that the susceptibility x 
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FIG. 8: (Color online) Bulk viscosity due to strangeness equilibration for the CFL-ZC" phase (solid lines), the CFL phase 
(dashed lines), and unpaired quark matter (dashed-dotted line). For all curves, the quark chemical potential is — 400 MeV, 
the compression frequency lo/{2tt) = lms~^, and the kaon mass m — 10 MeV. We have chosen different values of the kaon 
chemical potential /i. Positive values of 5m = m — fi correspond to the CFL phase (no kaon condensation) while negative values 
of Sm correspond to the CFL-ZsT" phase (kaon condensation). 



is discontinuous at Tc (and large for temperatures below and close to Tc) . Consequently, we expect the bulk viscosity 
to be discontinuous too (and small slightly below Tc). We show the numerical result for the bulk viscosity around 
the critical point in Fig. [T] The case of an almost continuous behavior as well as the case of a strongly discontinuous 
behavior is shown (for the latter we have chosen a frequency uj /(2tt) — 3 ms""'^ which is larger than any known compact 
star rotation rate). 



3. Parameter dependence and comparison to other quark phases 

Finally, in Fig. [SI we present the bulk viscosity for several values of the kaon chemical potential /i and compare 
the result with the bulk viscosity of unpaired quark matter [13]. As discussed in Sec. Ill Al the values of /i and the 
kaon mass m are poorly known at densities relevant for compact stars. We therefore vary the relevant parameters 
over ranges of values that are plausible for conditions in a compact star. We fix the kaon mass at m = 10 MeV and 
study six values of the effective kaon chemical potential /i; for three of them 5m = m — jjL la positive (corresponding 
to the CFL phase) and for the other three 5m is negative (corresponding to the GF\j-K^ phase). We have restricted 
our plot to temperatures appropriate to compact stars, i.e., we have not shown temperatures larger than 15 MeV. 

In the CFL phase, the energy gap for the kaon is 5m, so the thermal population of kaons, and hence the kaon-decay 
contribution to the bulk viscosity, will be very sensitive to the value of 5m, and will drop rapidly as exp(— 5m/T) 
for T <^ 5m. In the CFL-if" phase, on the other hand, there is always a massless Goldstone kaon, so the bulk 
viscosity due to kaon decay should be less sensitive to dm, and should not drop exponentially at low temperatures. 
These expectations are borne out in Fig. [5] Although the CFL-7^° phase, thanks to its Goldstone mode, has a higher 
bulk viscosity at very low temperatures, we see that in the range 10 keV < T < 10 MeV the CFL phase has a larger 
kaon-decay bulk viscosity than the CFL-K^ phase. This is because of the phase space available at the <-> H 
"resonance" (see the end of Sec. IIII Cp . There is another contribution to the bulk viscosity, from H ^ H + H, which 
starts to become comparable to the contribution from kaon decay at T < 1 MeV [l^] . 

The bulk viscosity of CFL/CFL-iiT" matter is only comparable to that of unpaired quark matter at relatively high 
temperatures, of order 10 MeV. At lower temperatures T < 5 MeV the bulk viscosity of CFL/CFL-i^" quark matter 
is several orders of magnitude smaller. This is the case not only in comparison to the unpaired phase but to all 
other color-superconducting phases with ungapped fermionic modes — these have bulk viscosities comparable to that 
of the unpaired phase, and even larger at high temperatures [2^. The only other color-superconductin g ph ase besides 
the CFL / CFL-i^° phase in which all fermions may be gapped is the color-spin locked phase [s^ [sol l69l. F70| . For a 
discussion of its bulk viscosity see Ref . [2l[ . 



17 



IV. SUMMARY AND CONCLUSIONS 

We have computed the bulk viscosity of kaon-condcnsed CFL quark matter (CFL-A'*^ phase). Kaon condensation 
affects the low-energy properties significantly and therefore has a significant effect on thermodynamic and transport 
properties of color-flavor locked quark matter. In particular, the CFL-isT" phase has a massless bosonic excitation 
associated with kaon condensation which is absent in the pure CFL phase. In both CFL- if ° and CFL phases, there 
is also a massless Goldstone mode H associated with superfluidity. For most of the thermodynamic properties, the 
effect of the additional Goldstone mode is important, but rather easy to predict. We have shown this for the specific 
heat, which acquires a contribution of the kaon mode which, as expected, has the same temperature dependence as 
the contribution of the superfluid mode. The prefactor of the former, however, is typically larger than that of the 
latter such that the kaon contribution in fact dominates the specific heat at low temperatures. 

The effect of the additional Goldstone mode on the bulk viscosity is more complicated. We have used the results 
of our earlier work [4^ which provides a self-consistent description of the CFL-i^T" phase for arbitrary temperatures. 
Using the resulting thermal kaon mass and excitation energy, we have computed the density and the susceptibility 
of kaons, both of which are needed for the bulk viscosity. Moreover, we have computed the rate of the processes 
^ H + H and + H ^ H, where we denote by if" both the neutral kaon in the CFL phase and the massless 
Goldstone mode arising upon kaon condensation in the CFL-if ° phase. These weak processes serve to re-establish 
chemical (flavor) equilibrium in response to an external volume oscillation, hence giving rise to a nonvanishing bulk 
viscosity. 

At very high temperatures, T > lOMcV, the difference in the kaon excitations in the CFL and CFL-if° phase 
is negligible. Consequently, in this case the kaon production (and annihilation) rate is almost identical for the two 
phases. At smaller (but not too small) temperatures, 10 keV ^ T < 10 MeV, the masslessness of the Goldstone mode 
in the CFL-if" phase suppresses this rate because of a smaller available phase space for the weak process. Since the 
timescale of the rates in both phases is smaller than the typical oscillation (and rotation) frequency in a compact 
star, this effect decreases the bulk viscosity of the CFL- if " phase compared to the CFL phase. Another effect is given 
through the different susceptibilities. In the condensed system, the susceptibility at low temperatures is much larger 
than that of the uncondensed system. This effect works in the same direction, further decreasing the bulk viscosity 
compared to the uncondensed system. For even smaller temperatures, T < 10 keV the phase space actually is larger 
in the CFL- if " phase and consequently the bulk viscosity is larger too. It is interesting to note that for the neutrino 
emissivity the effect of kaon condensation is quite different: neutrino emissivity in the CFL- if phase is larger than 
in the CFL phase for all temperatures T < 10 MeV [7l|- 

We now have a fairly complete understanding of bulk viscosity in color-flavor-locked phases of quark matter. 
The suppression of the bulk viscosity due to the absence of ungapped fermionic excitations was predicted in Ref. 
[l6| . Subsequent more careful calculations took into account the contribution of the superfluid [l^l and kaonic 24 1 
Goldstone modes. With the result of the present paper we have shown that the conclusion already drawn in Ref. [16 1 



is, for temperatures T < 1 MeV, not changed by the contribution of the Goldstone modes: color-flavor locked quark 
matter, even in the presence of kaon condensation, has a much lower bulk viscosity than all other known phases of 
dense quark matter and than nuclear matter. Only at large temperatures, and thus in very young neutron stars, can 
the contribution of the Goldstone modes render the bulk viscosity comparable to that of unpaired quark matter. 

We finally mention that besides the bulk (and shear) viscosity, other properties of color-superconducting quark 
matter also deserve attention. Its equation of state may be used to put constraints on the mass-radius relation 
of hybrid stars with a quark matter core and a hadronic mantle. These calculations rely on simple models whose 
parameters are poorly known in the strong-coupling region of interest. While NJL model calculations, mainly due to 
their relatively large predicted strange quark mass, tend to find no stable hybrid star with a CFL core [zl, [TJ], other 
parametrizations of the equation of state allow for hybrid stars with masses compatible with the observations [zl] . 
Other observables that may distinguish between certain phases of color-superconducting quark matter or between 
quark matter and nuclear matter are for instance the cooling curve of the star or glitches (sudden spin-ups). The 
corresponding transport properties of color superconductors have already been computed in the literature, see for 
instance Refs. [53,[6l,[7l,l7l] and Ref. [t^] for neutrino emissivity and shear modulus, respectively. It is an interesting 
and promising task for the future to extend these calculations and compare them with more and better astrophysical 
data in order to understand matter inside a compact star, and, ultimately, map out the phase diagram of cold and 
dense quark matter. 
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APPENDIX A: SPECIFIC HEAT IN THE 2PI FORMALISM 



Here we derive the expression for the specific licat in the 2PI formalism. We have to take into account that the 
condensate (j) as well as the full propagator S arc implicit functions of the temperature T (remember that S contains 
the self-consistent thermal mass M{T)). With the definition ((3T|) we thus have to compute 



cv = TV^P , 



under the contraint that at the stationary point 



V-—— — 
" dfdS ^ dTd^ ~^ df 



fs = 



^=0 
dS ' 



(Al) 

(A2a) 
(A2b) 



We find 

dT 



dT as dT 



dfs 



From the constraint (IA2D we obtain 



dT 



dT 



dSdU 
dT dS 



dcf) df^ 
dT d(h 



dT 



d^dfs 
dT dT 



d(f) df^ 
dT dT 



d^P 
gf2 



(A3) 



95 

df'df 



dfs 
dT ' 



dT 



( dfs_ dU\ 

dS dS 

dfs_ dU 

dcj) d<j) 



(A4) 



Inverting the matrix explicitly and inserting the result into Eq. (jASp shows that the expressions in parentheses vanish 
separately, and we are left with 



Cv 



TV 



dP 
'df 



(A5) 



We thus have to apply the differential operator T> on the explicit derivative of P with respect to T. From Eq. 
we find 



dP 
df 



(2^ 



1 



{[1 + /(efc)] ln[l + /(cfc)] - /(£,) ln/(efc)} + -(M^ - m 



2a(i 



2al) 



d]_ 

df 



(A6) 



At the stationary point, the second term in this expression vanishes. This leaves the first term as the result for the 
entropy density. In order to compute the second derivative, however, we have to keep the second term. By applying 
the differential operator V on this expression, the specific heat can be evaluated in a purely numerical way. For a 
further analytical evaluation we proceed as follows. In terms of the self-consistent quantities (j) and M we can write 
Cv as (since P only depends on the squares 0^ and Af ^ we may take the derivatives with respect to the squares) 



Cv ^T 



f d^P 



d^P dAP 



\dTd4P- dT dTdA'P dT 



d^P \ 
dT^ J 



(A7) 



The second derivatives of the pressure are straightforwardly obtained from Eq. (jA6p . The derivatives of (fP and AP 
with respect to the temperature depend on the stationarity equations. Below the critical temperature, one obtains 
the derivatives from Eqs. (|14a[) and (jl4bp . Above the critical temperature, Eq. (|14a[) is automatically fulfilled and the 
only equation is Eq. (jl4bp with (j> set to zero. Consequently, the derivatives assume different functional forms below 
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and above the critical temperature. We find 

_ 

dT 



dT 



dl 



2a- 



dl 



dT \ dA'P 



dl 



dT \ dhP 



-2a— 2a 



for T <Tc 



for T > Tc 



1 for T < Te , 



(A8a) 



(A8b) 



„ 5/ , ^ 
2a— 2a 



dl 



dT \ dhP 



1 for T > Te , 



Now we can insert these expressions and the second derivatives of the pressure obtained from Eq. (|A6p into Eq. (|A7I 
and evaluate the result at the stationary point, i.e., we use Eqs. p^ . The result is 



1 



T dl dhP 

2 J (27r)3 r2 cosh ^ - 1 ^ 2"9r dT 



(A9) 



With Eq. (jA8bp one then obtains the results below and above the critical temperature. The discontinuity bXT ~Tc 
then is 



Acy = cv{T^ T+) - cy(T ^ T-) = 2aT 



dl_ 

dT 



1 



2 ■ 



(AlO) 



APPENDIX B: BULK VISCOSITY WITH QUARK NUMBER EFFECTS 



In this appendix we derive the bulk viscosity taking into account the oscillations in the quark chemical potential 
as given in Eq. (|43bp . Instead of the simplified single differential equation this yields the two coupled equations 



dUq^djiq^ ^ 
d^lq dt 
dn d^q 

dflq dt 

Wc insert the form of the volume oscillation 



Im 5fiq 
Im 5fi 



duq dfi 
d^i dt 



duq dV 
Wdt 



(Bla) 
(Bib) 

and Eqs. into these differential equations and find the solution 

SVo H'^^-^nq) dUq 



-UJ- 



Vo (detJ)2w2 + (|^A)2 9^ ' 
5V^ dnq 



(det J)2w2 + (|1^A)2 duq ' 
where J denotes the Jacobian of the two- valued function [nq(/ig, /i), n{^q, ^)], 



(B2a) 
(B2b) 



J 



/ dUq^ dUq^ \ 

d^q dfi 

dn dn 

V d^iq dfi J 



(B3) 



Inserting Eqs. (|B2p into the definition of the bulk viscosity yields 

c 



d^j, d^lq 
dn dn 



d^ dfi 
n— h n, 



dn ' dn.. 



on 



(B4) 
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This is the general form of the bulk viscosity where the derivatives of the chemical potentials with respect to the 
densities are obtained by inverting the Jacobian (|B3p . We now assume that the quark chemical potential enters the 
densities only through the kaon chemical potential, see Eq. ([5]). This means that we neglect the dependence on 
through and a in Eq. Defining the small dimensionless quantity 
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(B5) 



we then simply have Uq = —nrj and 



J 



XV 



2n 

— -n 



-XV 



-XV 
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dflq_ 


djlq \ 


( 




dUq 


dn 


_ 1 2nri 






dix 




V 


dUq 


dn ) 


V 2n 



2n 



TT X 



(B6) 



Consequently, Eq. (|B4p becomes 



[{x-'+^v)xr 



(B7) 



Neglecting the term proportional to rj yields the result (|46p that we use in the main part of the paper. 



APPENDIX C: NUMBER SUSCEPTIBILITY IN THE 2PI FORMALISM 



In this appendix we derive an expression for the kaon number susceptibility x. Since x is given by the second 
derivative of the pressure (|18p with respect to the chemical potential, we can use the same formalism as in Appendix 
1X1 where we have computed the second derivative with respect to the temperature. Consequently, with the same 
arguments as in Appendix El (cf. Eq. (|A7p ) 



d^P dAP d^P 



We determine the derivatives of 0^ and with respect to /i from the stationarity equations ((H]). They are 



(CI) 



r 2/i 

Ofj,a 

a a 



2 a {d^I + AfidMl) + Id^a 



1 + 2adMl 







for T < , 
for T>Tc, 



(C2a) 



An- 2 



1 + 2adMl 
1 - 2adMl 



for T <Ta 



for T>Ta 



(C2b) 



where we have abbreviated 9^ = d/dfi, Om = d/dM^. Note that the effective coupling constant a depends on the 
kaon chemical potential ^, see Eq. Next, we compute the second derivatives of the pressure appearing in Eq. 
(|C1[) with the help of Eq. (|18p . Putting everything together we obtain the susceptibility. 



Xo + Xi + 



2T J (27r)3 cosh ^ - 1 



2T J (27r)3 cosh ^ - 1 



for T < Tc 



for T>Tr 



(C3a) 



where the zero-temperature part is 



Xo 



a 



a4 

4 



(C4) 
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and where 



/a 



1 + 2adMl V a 1 + 2aaM/ 

+ + (C5a) 

, _ adf,I + Id^a 



APPENDIX D: LOW-TEMPERATURE RESULT FOR THE KAON RATE 



In this appendix we derive Eq. (|62p . the low-tcmpcrature expression for the kaon rate in the presence of a kaon 
condensate. To this end, we first compute the low-temperature result for the imaginary part of the retarded H 
self-energy. Below the critical temperature we have Cp < vp for all p (see Eq. (HH)) and thus we have 

Imn(ep,p) = n-(ep,p) = ^g^^^4^7 pj(^,^) dkg\ep,p,k)fivk)[l + f{vk - e^)] . (Dl) 

For small temperatures only small momenta p < T are relevant, and we can use the small- momentum approximation of 
the excitation energy given in Eq. (|20p . Moreover, since we consider temperatures well below the critical temperature, 
the quadratic part of the kaon dispersion is negligible. Therefore, ~ (3vp, with the abbreviation 



f3^f3{m,n)^J f\ (D2) 

Here we have used the zero-tcmperature result for the self-consistent mass M . With the new integration variable 
X — vk/T we have approximately 

Imn(ep, p) c dx - 1)^;V - 3(1 - r^)Tx{f3vp - Tx)]^ . (D3) 

Note that p and T are of the same order. Therefore we cannot simply drop the terms in the square brackets that are 
of higher order in T . We rather have to keep all the terms. The integration can be done exactly and we obtain 

Imn(ep, p) ~ ""flf e-^^ [{(i^ - 1) V + 8^3 {(3^ ~ 1) p^T + 48 {[3^ + 1) p^T^ + 576V^pT' + 3456 T"] , (D4) 
21b/ig/ (Epj 

where, for the sake of brevity, we have used v = 1/^/3- We can now insert the result into the rate (|60p . In the 
denominator of the integrand we neglect the term Imll^. No further approximation is then required to obtain the 
final result. We make use of 

, x"e-^ r(n-t- 1) / 1 + , s 



sinh/3a; 2"/3"+i V 2/3 
where C is the generalized zeta function, and where we need the cases rt = 2, . . . , 6. Then we arrive at 

with the dimensionless function 

Inserting the expression for f] from Eq. (|D2p into Eq. (|D6|) yields the final result given in the main text, see Eq. ((^ . 
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